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Abstract — The goal of this paper is to provide conditions 
under which a quantum stochastic differential equation (QSDE) 
preserves the commutation and anticommutation relations of 
the SU(n) algebra, and thus describes the evolution of an open 
n-level quantum system. One of the challenges in the approach 
lies in the handling of the so-called anomaly coefficients 
of SU(n). Then, it is shown that the physical realizability 
conditions recently developed by the authors for open n-level 
quantum systems also imply preservation of commutation and 
anticommutation relations. 

I. Introduction 

The property of physical realizability of open quantum 
systems has attracted considerable interest in recent years 
with its main motivation being its role in quantum coherent 
control [1], [2], [11], [15]. In simple words, this property 
systematically characterizes the quantum nature of a control 
system from a state space point of view, which is relevant in 
many engineering areas, and is closely related to algebraic 
properties of the QSDE coefficients. The case of open quan- 
tum linear systems is well understood [8]. More recently, 
the physical realizability property was introduced for ?i-level 
quantum systems [4]. This was the product of a nontrivial 
extension of the two-level case presented in [3]. From the 
state space point of view, physical realizability has to be 
complemented by the preservation of commutation relations 
since the latter is an exhibitor of the quantum behavior of 
the system. In this regard, the authors provided conditions 
for the preservation of commutation relations in the case of 
QSDE's evolving in SU(2) [5]. However, in order to extend 
the formalism to the general n-level case it is necessary to 
deal with the so-called anomaly coefficients of SU (n), which 
forms the completely symmetric tensor dijk that appears 
out of the anticommutation relation of the generators of 
SU{n). These generators form a complete orthonormal basis 
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spanning the set of n-dimensional complex matrices and are 
known as the generalized Gell-Man matrices [10], [12], [13]. 
The commutation relations of these generators also generate 
the completely antisymmetric tensor /y^. In the case of 
QSDEs describing two-level quantum systems evolving in 
SU(2), it suffices to consider preserving the commutation 
relations of SU(2) since the anomaly tensor d is zero for 
all indices. On the other hand, the situation is nontrivial for 
n > 3 since one now requires extra machinery on the tensor 
d that in some cases depend on the value of n. In other words, 
the QSDE for an n-level quantum system must preserve both 
commutation and anticommutation relations for all times. 

The approach followed in this manuscript is that of open 
quantum systems defined in terms of a triple (S,L,H), 
where S is an operator, known as the scattering matrix, 
describing the interaction of the environment fields among 
themselves, L is a vector of coupling operators expressing 
the interaction of the environment fields with the system 
variables and H is a Hamiltonian operator [6], [9]. This 
description implicitly takes in consideration the quantum 
nature of the evolution equations. For example, the (S, L, H) 
description will preserve, in time, the canonical commutation 
relations of the system under consideration. The physical 
realizability conditions in [4] do not, in general, assure the 
preservation of commutation and anticommutation relations. 
It is thus the main goal of this paper to provide conditions 
for the preservation of commutation and anticommutation 
relations for a state space model evolving in SU(n) inde- 
pendently of the physical realizability conditions. Then, as 
a second result it is shown that physical realizability does 
imply the preservation of commutation and anticommutation 
relations. 

The paper is organized as follows. Section HI1 summarizes 
the necessary tools obtained from the algebra of SU(n). In 
Section [III] the basic preliminaries of open ?i-level quantum 
systems as well as the property of physical realizability for 
such systems are introduced. This is followed by Section 
IIVI in which conditions for the preservation of commutation 
and anticommutation relations are provided. In addition, 
it is also shown that those conditions are implied by the 
physical realizability conditions. Finally, Section [Vl gives the 
conclusions. 

II. The algebra of SU(n) 

In what follows the necessary tools regarding the algebra 
of the special unitary group SU(n) will be provided; see 
[10], [12], [13] for more details. Basically, this group is 
formed by all n x n complex matrices that are Hermitian 
and have zero trace. Consider the set of elementary vectors 
spanning C", namely, {ei, ■ • • , e„}. Define P^i £ C" x ™ as 



Pk.i = Gkef, where k,l = l,...,n. A standard way of 
constructing a complete orthonormal basis for SU (n) is 

Ujk = Pj.k + Pk,j> 
Vjk = i (Pj,k - Pk,j) i 



Wl 



for 1 < j < k < n, 1 < I < n - 1. Note that the 
identity matrix I must be included in the basis in order to 
form a complete set. The identity, (n 2 — n)/2 symmetric 
matrices Uj k , the (n 2 — n)/2 antisymmetric matrices Vj k and 
the n — 1 mutually commutative matrices wi together form 
the generators of SU(n). These generators are known as 
the generalized Gell-Mann matrices. Without any particular 
order, the generators are relabeled {/, Ai, . . . , X s }, where 
s = n 2 — 1. Here the orthonormality condition these matrices 
satisfy is Tr(XiXj) = 2<5y, where <5y denotes the Kronecker 
delta. Their commutation and anticommutation relations are 

8 

[A,, Xj] = 2i fjjkXk, 

k=l 

{Ai,Aj} = —Sij + 2 ^ djjkXk- 

n k=l 

Thus, the product X{Xj can be easily computed as 

AiAj = -([X i ,X j ] + {X i ,X j }) 



% + / , + dijk) X k 



(1) 



where the real completely antisymmetric tensor and the 
real completely symmetric tensor d^k are called the structure 
constants of SU(n). The tensors f^k and dijk satisfy 
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film^mjk fjlm^imk fklmdijm 0; 



k=l 



f ilk f mjk — \^im^lj ^ij^lm) 



m.k—1 



+ ^ (dimkdljk 
k=l 

fimkfjmk — TlSij . 



^ijk^lmk) i 



(2a) 
(2b) 



(2c) 



(2d) 



Define F ? , A G 



i G {1, . . . , s}, such that their (j, k) 



component is (Fi) jk = f ijk and (A)jfc = <%fe, respectively. 
In particular, the set {— iF\, . . . , — iF s } is the adjoint repre- 
sentation of SU(n). In [10], [12], identities (|2a|)-(|2c| were 
employed to obtain the following useful relationships 



[Fi,Fj] = - fi 3 kFk 

k 
s 

[F i> D j } = -J2fmD k 

k 

s 

F l D J + FjDi = ^d ijk F k 



(3a) 



(3b) 



(3c) 



DiFj + DjFi = ^d ijk F k 

k 
s 

{DiDj -FjFi) ml = Y^d ijk (D k ) 



(3d) 



H (SijSml — 5i m Sjl) . 



(3e) 



Definition 1: Let /3 G C s . The linear mappings 9 , + 
s — >• C sxs are defined as 

( P T F? 

e-(/3) = (F 1 T /3,-.-,^ T /3) = 



e+C8) = (Dfp,.-- ,i%p) 



Observe that the nature of the / and rf-tensors make 0~(/3) 
and + (/3) be antisymmetric and symmetric, respectively. 
When j3 is an s-dimensional row vector then it will be 
understood hereafter that 9~(P) = e~(P T ) and 9+(/3) = 
G + (/3 T ). Consider now the stacking operator vec : C mXTl — > 
<C mn whose action on a matrix creates a column vector by 
stacking its columns below one another. With the help of 
vec, the matrices 9~(/3) and + (/3) can be reorganized so 
that 

/ eros) 

vec(6-(/3)) = : 

and 

vcc(Q+(P)) = : 

where /3 G C 5 , 0r(/3) = Ff/3, F 
Of (/3) = A/? and D = (D x , ■ ■ ■ ,D S ) T . From <T2dJ, F 
satisfies 

F T F 



F0, 

Dp, 
(Fi, 



nl. 



(4) 



The properties of 6 and + are summarized in the next 
lemma. 



Lemma 1: Let /3, 7 G 
satisfy 

e-(/?) 7 = -e-(7)/3, 
e+(/3) 7 = e+( 7 )/3, 

e- (e-(/?) 7 ) = [e-(/?),e-( 7 )] 



The mappings 9 and + 



e- (e+(/3) 7 ) 
e+ (e-(/?) 7 ) 
e+ (e+(/3) 7 ) 



e-(/3)e + ( 7 ) + e-( 7 )e + (/3), 

[e+C8),e-( 7 )] = [e-( ) 8),e+(7)] 
e + (/3)e+( 7 )-e-( 7 )e-(/3) 

- - (/3V - /? 7 T ) ■ 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 
(5f) 

(5g) 



Proof: The proofs of identities (l5all-(l5fTl where provided by 



the authors in [4]. Only (5gi requires a proof. The left-hand- 



side of ( pg| ) is decomposed as 

e+ (e + (/?) 7 ) 

- U { ,,TDl1 

I 



I ? ^ s 



( P T D ll 



k=l 



k=l 



k=l 



\ fc=i 

By d3eb . the (i, j) component of this matrix is 

(e+ (e + (/3) 7 ))„. 



\k=l 



:/3 T (I), I), I J-,)-. 



m,l — l m,l — l 

2 



= r D i D i 7 - !3 T F 3 F tl - - (/3 T 7 <% - ft 7i ) 

= (e+(^)e+(7)-e-(7)e-(^)).. 
-!(^V-^7 r L. 



Some additional identities regarding matrices F and D 
with respect to the Kronecker product are given next. They 
will become useful when proving the main results of the 
paper in Section [IV] Define the tensor permutation matrix 
1® G R s xs . That is, a symmetric block matrix 1® = 
{ljt}f j=i such that it satisfies 1®(A ® 5)1® = (-B ® A) 
for any A, B e C sxs . 

Lemma 2: Let F = F 2 , • • ■ , F S ) T and F = 

(D U D 2> --- ,D S ) T , and A, B E R sxs . Then 



i. F = -1 9 F, 

a. d = i®d, 

Hi. F T (A <g> F)F = F T (5 (8 A)F, 
w. D T (A ® 5)1) = D T (B ® A)D. 

Proof: Identity (0 is proved by multiplying directly 1® and 
F. Therefore, 

/In ••• l.i \ / Fi N 

I / 

V lis ••• la. / V/s 




\ \( F s)sJ J 

Observe that (F^j = (fiji,--- ,fij s ) = ~(F)j since f ijk 
is an antisymmetric tensor. Using this fact in the previous 
expression gives 



1«F 



FT 



F 



= -F 



An analogous procedure is used to show identity (Jm|), but in 
this case the permutation of indices in dijk does not produce 
a negative sign due to the fact that dijk is a completely 
symmetric tensor. For (|m[) . one will employ (0) as follows 

F T (A ® S)F = (-1^F) T (A g> S)(-1 8 F) 

= F T (1^(A®B)1 )F 

= F T (B (gi A)F. 



Finally, identity is proved similarly but using (|m|) 

instead. ■ 



III. Physical Realizability of Open ti-Level 
Quantum Systems 

Consider the separable Hilbert space f), T(i5) the set of 
operators in fj, and c I(S)) nxrn the set ofnxm dimensional 
arrays of operators in T(-fj). The commutator of x and y in 

1(Sj) is [i, y]= xy — xy. Let i G T(io)" 1 and y G T(^)" 2 , 



T] A 



(yx 7 ) 2 G T(io)" lX " 2 . The adjoint of 




then [x, y 
x is denoted by x 1 



and * denotes the operator adjoint. For the case of complex 
vectors and matrices, * denotes the complex conjugate while 
t denotes the conjugate transpose. 

The Heisenberg evolution of a; £ T(-fj) s interacting with 
n w quadrature Boson quantum fields W\ and W% in the 
parametrization (S,L,H) is given by 

dx = C{x) dt+- (jx, L T ] - [x, L f ]) dWi 



% -([x,L T ] + [x,L^])dW 2 



(6) 



where C(X) is the Lindblad operator defined as 
C{x) 

= -i[x,H] + \((L^ [x,L T ] T ) T + [L#,x T f 
The output field in its quadrature form is 

'd¥i\ ( L + L* \ fdWi 
: {i(L#-L)J M+ \dW 2 



dYi 



dt. (8) 



Here the operator S is assumed to be the identity operator 
(S = I), and the Ito table ( [7]) for the fields W 1 and W 2 is 

(*) ( « m ) . ( x "z 

The interest in this paper is in systems evolving with 
respect to the special unitary group SU(n). Consider the 
Hilbert space for these systems to be = C n . It is standard 
to associate a vector /3 = - ,(3 S ) G C s with the 

vector of operators j3 G 1(f)) s by simply considering the 
identity operator in each component of the vector, i.e., j3 = 
(/?!/,••■ ,/3 s I). In a similar manner, any complex matrix 
is associated to a matrix of operators by considering the 
identity operator I in each component. The identity operator 
I is usually suppressed. Therefore, abusing of the notation 
slightly, the mappings 0~ and 9 + are allowed to act on a 
vector of operators. In this manner, it is valid to multiply 
complex matrices and operator matrices. This convention 
allows to express the commutation and anticommutation 
relations of the system variables x as 

[x,x T ] = 2ie-(x), (9a) 
4 



{x,x T } = -7 + 26+(x). 



(9b) 



The vector of system variables for (|6) evolving in SU(n) is 
then 



where Ai, . . . , A s are self-adjoint and spanned by the gen- 
eralized Gell-Mann matrices. They are usually called spin 
operators. The initial value of the system variables can 
be set to x(0) = (Ai,...,A s ) with Ai,...,A s being the 
generators of SU (n) introduced in Section |TjJ Due to the 
product relation (Q3 any polynomial of spin operators can be 
written as a linear combination of generalized Gell-Mann 
matrices. Therefore, assuming linearity captures a large 
class of Hamiltonian and coupling operators. This allows 
to assume, without losing generality, that the Hamiltonian 
is H — ax with a G R s , and the multiplicative coupling 
operator is of the form L = Ax with A G C n ™ xs . 

In general, the evolution of x in quadrature form falls into 
a class of bilinear QSDEs expressed as 

dx = Aodt + Ax dt 



(B n x,--- ,B lnm x, B 21 x, ■■■ ,B 2nm x) f^p*^ > 



(10) 



dYy 

d% 



(11) 



where A G W, A, B lk 4 B lk +B 2k ,B 2k 4 i(B 2k -B lk ) G 
R sxs and C 1 ,C 2 G r» xs , k = l,...,n w . It is assumed 
that all matrices in ( TTOb and (fTTl i are real due to the fact that 
the class of quantum systems considered in this paper are 
in quadrature form. The next results in physical realizability 
were given in [4] and are employed later in Section II VI 

Definition 2: [4] A system described by equations (TlOb 
and (fTTb is said to be physically realizable if there exist 7-L 
and L such that ( TTOb can be written as in ([6j. 



The explicit form of matrices Aq, A, Bi k ,B 2k , C\ and C 2 
in terms of the Hamiltonian and coupling operator is given 
next. 

Theorem 1: [4] Let % = ax, with a T G K s , and L = 
Ax, with A G C n ™ xs . Then 

Ai " 



k=l 

yi in 

A= -2Q-{a)+Y^(R k -iQ k ), (12b) 



fc=i 

B lk = Q- (i(A* - A fc ) 

B 2k = -Q-(A k +Af), 
d =A + A# 
C 2 = i (A* — A) , 



(12c) 

(12d) 
(12e) 
(12f) 



where 







/AA 


i. A = 








a. Bi k = 



R k 4 e-(A fe )e-(A#) + e-(Af )e-(A fc ), 
Q k 4 e-(A fc )e+(A#) - e-(A# )e+(A fc ). 

The next theorem present conditions for physical realiz- 
ability in terms of (Aq, A, B iy C) in (|T0j and (fTTb . 

Theorem 2: [4] System ( fTOt with output equation (fTTb is 
physically realizable if and only if 

i «™ 

A) = - y"(»Sifc + B 2k ) ((Ci) fc +i(C 2 ) fc ) T , 

fe=l 



^e+(A ), 



m. S 2 * = e-((<7i) Jk ), 
i«. A + A T + B ik B lk 

i,k=l 



where (Ci) k indicates the fc-th row of Ci. In which case, the 
coupling matrix can be identified to be 

A = ±(C 1 +iC 2 ), 
and a, defining the system Hamiltonian, is 

a = -vcc [A T A + - Y, {[B 2 k, e + ((C 2 ) fe )] 
n \ fc=i 



[ J Bi fe ,e + ((c 1 ) fe )] 



(13) 



IV. Preservation of Commutation and 
Anti-Commutation Relations 



In this section, the main results of the paper are provided. 
It is shown first under what conditions a state space model, 
described by the QSDE ([Toi l, preserves the commutation and 
anticommutation relations of SU(n). The employed proce- 
dure for this task is independent of the physical realizability 
conditions given in Theorem [2] Also, recall from Section Hill 
that G G C sxs can be always regarded as G G C Z(S)) SXS since 
the identity operator in Z(Sj) can be considered attached to 
each component of G. The following Lemma plays a key 
role in the forthcoming theorems. 

Lemma 3: Let G G 1(Sj) sxs be antisymmetric and x 
be a vector comprised by elements of a complete linearly 



independent set of operators spanning 'X(io). If G satisfies which means that Gy = — ^ J2k=i fijk^(FkG). Then, 

Ge-(x) + e-(x)G T -@-(Gx)=0, (14) _1A 
then there exist g£%(Sj) such that ij n jk ' k ' 



G = e-(g), (15) /Tr(F x G) N 

where the unique g is given by = . . . ; fj is ) 



5 



A 1 



[MFiG)\ \Tr(F s G) / 

/Tr(F!G) N 



V Tl '(^G) / 



1 rp 

"7-1 



/Tr(FxG) N 
^Tr(F s G) / 



(16) 

v*m/ =-^)ii,---,(^ s ) 

Conversely, if < fT3T > holds, then dT4l > is true for any x G 1(Sj) s . 

Proof: Assume (O holds. By applying vec to (O and Therefore, 
recailing that F = (Fi, • ■ • , F S ) T , it follows that 

vcc (Ge~(x) + Q-{x)G T - Q-(Gx)) G 3 = ~Fj 
= {I® G)Fx + {G® I)Fx ~ FGx. 

Since by assumption the components of x are linearly B Y defining g as in (16), one obtains that 

independent, it can be concluded that Q = (F^g, ■ ■ ■ ,Fjg) = 0~(g). 

(J (g> G)F + (G <g> J)F - FG = 0. (17) Now, assuming that there exist a 5 such that G = 8~ (.9), 

The application of identity (Jm) in Lemma [2] and © after il then follows directly from (gd) that the left-hand side of 

multiplying on the left by F T gives £0} satisfies 

« _ e-(g)e~(x)-e-(x)&~(g) T -e-(e-( s )x) = o 

fc fe + n ■ £ or an afbitrary x £ which concludes the proof. ■ 

The equation for the component of G is then In order to be considered a quantum system, the system 

s variables of ( [Tol l must preserve the product rule (H) for 

2 fkjrfkuGir — nGij = 0. (18) all times. This amounts to preserve the commutation and 

k,r,i=i anticommutation relations of SU(n). Thus, the conditions 

From @ and noting that Tr(AB) = ££ J=0 A rl B lr for any that CE1 must satisfy are 

A,Be IT* 3 , O) is rewritten as d[x, x T ] - 2iQ~ {dx) = 0, (20a) 

n ^ d{x,z T }-2e+(dx) =0. (20b) 

U — ^ y Jkjr Jkil^lr — nt^rij 

k l_ x Note by the linearity of 9~ and 6+ that 

Q~{dx)= Q-{A Q )dt + Q-{Ax)dt 

= - 2 Gw YUjikhkr + frikfijk) - nGij + Q-(B lX )dWi + 0-(B 2 x)dW 2 , 

r ' l= ? Q + (dx) = e + (A a )dt + Q + {Ax)dt 

= -2 J2 fjikhkr{G T ) r i +e + {B 1 x)dw 1 + e+{B 2 x)dw 2 . 

k - r ' l=1 A condition for system ( TTOt to satisfy ( I20al i and (I20bb is 

, r> r< given in the next theorem. 

- l 2^JHh 2^ JrlkGir - nGij Theorem 3: Let x(0) <E T(i5) s be comprised by the 

~* generators of SU(n). Then, system ( TTOb implies 

= -2 J] fjikhkr{G T ) rl [x(t),x(t)]=2ie-(x(t)), 

{x(t),x(t)} = ^I + 2Q+(x(t)) 

- 2 E /yfeTr(F fe G) - nG i3 - . (19) f or a ll t > 0, if and only if 

ai . fe = 1 B< = _ (6<) (21a) 

Also ((18) implies that v ' y ' 



k,rj—l fc,r,Z— 1 



2 £ /j,*/*^),, = 2 £ fati-foX-Grl) =nGij. X>**& = 2 6 ~ (Ao) (21b) 

, n w 

A = Q-(a)--Y {B lk Bl k + B 2k B^ k ) 
Substituting this on ( fT9l l gives k=1 

= -ndj ~ 2 £ /^Tr(F fc G) - nG l3 , + ~ 2 E - Si*e + (M) . (21c) 

* * j i 



k=i fc=1 



where bik and a are s-dimensional vectors as in (fT6] l for 
i = 1 , . . . , s and k = 1 , . . . , n w . 

Proof: Without loss of generality, one can consider (TlOb 
interacting with only one quadrature field. That is, 

dx = A dt + Ax dt + BiX dW\ + B 2 x dW 2 . 

Using the Ito table ©, it follows that 

d(xx T ) = (dx)x T + x(dx) T + (dx)(dx) T 
= (A Q x T + xA^) dt + (Axx T + xx T A T ) dt 

+ (B lX x T + xx T B'[) dWi + (B 2 xx T + xx T Bl) dW 2 
+ Bixx 1 'Bj 'dWidWi + Bixx T BldW 1 dW 2 
+ B 2 xx T BfdW 2 dW 1 + B 2 xx T B 2 r dW 2 dW 2 
= (A x T + xAq) dt + (Axx T + xx T A T ) dt 

+ {B r xx T + xx T Bl) dWi + (B 2 xx T + xx T B^) dW 2 
+ B r xx T Bfdt + iB 1 xx T B 2 dt 
- iB 2 xx T B\dt + B 2 xx T B 2 dt. 
Similarly, 

(d(xx T )f 

= (A x T + xAl) dt + (A(xx T ) T + (xx T ) T A T ) dt 
+ (B 1 (xx T ) T + {xx T ) T B^)dW l 
+ (B 2 (xx T f + {xx T ) T Bl) dW 2 
+ B 1 (xx T ) T Bldt+iB 1 {xx T ) T B^dt 
- iB 2 {xx T ) T B'[dt + B 2 {xx T ) T B 2 T dt. 
Computing d [x, x T ] then gives 



d{[x,x T }) = d{xx T ) - (d(xx T )) 
2 



2i ( -BiB? - -B 2 Bl I dt 



n 



ii 



+ 2i i^AQ- (x) + 8" (x)A T j dt 

+ 2i^BiQ-(x)Bf +B 2 Q~(x)Bj^j dt 

+ 2i(^B 1 e + (x)B^ - B 2 Q + {x)B^ dt 

+ 2i (Bi@-(x) + e~(x)B'[) dWi 
+ 2i(B 2 Q-(x) + Q-(x)B 2 r )dW 2 . (22) 
Replacing ( f25b into J20at amounts to 

2i (B 2 Bf - B x Bl - Q-(Ao) + AG~(x) + Q-{x)A T 

+B 1 e-(x)Bf + B 2 Q-{x)Bl - ®-{Ax)) dt 

+ 2i {BxQ-{x) + Q-(x)B'[ - e-(B lX )) dWi 

+ 2i (B 2 @-(x) + Q-(x)B% - Q-(B 2 x)) dW 2 = 0. (23) 

From [14, Proposition 27.3], one can also equate the in- 
tegrands in (|23l to zero. Also, recall that that x(0) is 
represented by the complete orthonormal set constituted by 
the the generalized Gell-Mann matrices. This implies that 
any linear combination J2t=o a i x i(®) unless a.; = 
for all i and G C. In addition, no linear combination of 
generalized Gell-Mann matrices generates the identity. So, 
given that x(0) ^ 0, any equation involving the system 
variables of the form Ax = b (A G C sxs ) implies A and 



b must be identically 0. Thus, the equations to be satisfied 
for preservation of commutation relations are 

B t Q-(x) + e~(x)B'[ -e-(B-Lx) = (24a) 



B 2 Q-(x) + Q-{x)BT - e~{B 2 x) 
-B X B T 2 - -B 2 B{ - Q-{A ) = 



(24b) 
(24c) 

AQ-(x) + Q-{x)A T + B 1 Q-{x)Bl + B 2 e~(x)Bj 

+ B 1 Q + {x)Bl - B 2 e + (x)Bf - Q-(Ax) = 0. (24d) 

In the same vein for the anticommutator, one has that 
d ({x, x T }) is given by 

d({x,x T }) = d(xx T ) + (d(xx T )) T 

= - (A + A T + B x Bl + B 2 Bl) dt 
n 



2 [A x 1 -VxA^ dt 



+ 2 [AQ + {x) + <d+(x)A 1 



dt 



2 B l Q + {x)Bl + B 2 + (x)B 



dt 



( BiQ~ (x)B 2 - B 2 0-(x)Bf ) dt 
'2 2 



+ 2 l-Bi + -B{ + B^+ix) + Q+(x)Bi dWi 
\n n J 

+ 2 \-B 2 + -Bl + B 2 0+(x) + e+(x)Bj) dW 2 . (25) 
\n n ) 

Using the same argument employed for the commutation 
relation of x, the integrands in ( 1231 are equal to zero. Thus, 
for preservation of anticommutation relations one has that 

B^+ix) + Q + (x)B^ -Q + (B lX ) =0 (26a) 

B 2 e + (x) + e + (x)B 2 r -0 + (B 2 x) =0 (26b) 

A + A T + B 1 Bf + B 2 BT -^Q + (A ) = (26c) 

B l + Bl = B 2 + Bl = (26d) 
A x T + x Al + A6 + (x) + e+(x)A T - Q + (Ax) 

+ B 1 e+(x)Bf + B 2 Q+(x)Bl 

- B 1 Q-(x)Bl + B 2 e-(x)Bf = 0. (26e) 

Assume that d24at -( l24dt and (I26at - (l26et hold. Condition 
( |21bt appear explicitly in ( |24bt . From ( |24at , ( |24b) and ( |26dt , 
condition (121 al l is obtained by direct application of Lemma|3] 
i.e., Bi = 0~(bi) holds, where bi is given in ([TBI for 
i G {1,2}. In order to show ( |21c) , A need to be expressed 
in a more convenient form. This is achieved by using (l5dl i 
and ( |24bl . That is, 

Q-(A ) = - (e-(b 2 )Q-(h) - e-(b 2 )Q-(h)) 

n 

= -e- (Q-(b 2 )h) . 

Due to the linearity of 6~, Ao is then uniquely determined 
by A = ^Q~(b 2 )b 1 . Computing Q + (A Q ) in terms of b 1 
and b 2 using < T5fb gives 

e+(^ )=-e + (e-(6 2 )6i) 

n 



= - (e + (& 2 )e-(&!) - e-(6i)e+(6 2 ) 

n 

+ e-(b 2 )Q+(h) - e+(6i)0-(ba)) . (27) 
Using ( I26cb and d27l i. one obtains the following expression 

e-(6x)e-(&x) + e-(6 2 )e-(& 2 ) 

= A - \ (e-(6a)e+(6i) - 0-(6i)e+(6a)) 

+ ^-i(e-(b 2 )e + (6i)-e-(bi)e+(b 2 ))^ . (28) 

Now let P 4 A-i (e-(6a)e + (6i) - 0-(&i)0+(6 2 )). One 
can trivially decompose P in its symmetric and antisymmet- 
ric parts. That is, 

P=\{P + P T ) + \{P-P T )- 
Define P = \ (P - P T ). Note from that 

P + P T = e-(6 1 )e-(6i) + e-(6 2 )e-(6 2 ). 

Hence, P can be written as 

p = a - 1 (e-(6i)e-(6 x ) + e-(6 2 )e-(6 2 )) 
-i(e-(& 2 )0 + (6i)-e-(6i)e+(6 2 )) 

It is only left to show that P can be written in terms of 0~, 
The formula for P allows to calculate 0~(Px) as 

0"(Px) 

= Q~(Ax) - i 0- (0-(6 1 )0-(6 1 )o; + 0-(& 2 )0-(& 2 )x 

+ 0-(6 2 )0 + (6i)x - e-(6i)e+(62)x) 
= Q-(Ax) 

- \ (e~ (0-(6i)0-(&i)x) + 6- (6-(6 2 )8-(6 2 )x)) 

- ^ (e~ (e-(6 2 )9+(6 1 )x) - 0" (0-(6i)0 + (6 2 )a;)) 

= Qr {Ax) + e-(6i)0-(a:)e-(6i) + 0-(& 2 )0-(x)0-(& 2 ) 

- i (0-(6 x )0-(6i)0-(x) + e-(6a)e-(6a)9-(x) 

+ e-(x)e-(6i)e-(6i) + 0-(x)0-(6 2 )©-(6 2 ) 
+ 0-(6 2 )0 + (6i)0-(x) + e-(b 2 )e+(x)e-(b 1 ) 

- e-(b l )e + (x)Q-{b 2 ) - e-(x)Q + (h)e-(b 2 ) 

- e-(6 1 )e+(6 2 )e-(x) - ©-(6i)@ + (x)©-(6 2 ) 
+ e-(6 2 )e+( a; )e-(6i) + e-(a ; )e + (& 2 )e-(6i)) . 

In the previous calculation, the antisymmetry of 0~ was 
used as follows 

0- (0+(& 2 )x) = 0-(& 2 )0+(x) + 0-(x)0+(& 2 ) 

= - (0-(& 2 )0+(x) + Q~ (x)0+ (b 2 )f 

= e + {x)G-(b 2 ) + G+{b 2 )e~{x). 

Also, 

P0"(x) + Q-{x)P T = AQ-(x) + G~(x)A T 

- 1 (e-(6i)e-(6i)e-(x) + e-(6a)e-(6a)e-(x) 

+ 0-(x)0-(6 1 )0-(6 1 ) + e-(x)e-(6a)9-(6a) 
+ 0-(6 2 )0+(6 1 )0-(x) - e-(x)e + (6i)e-(6a) 



- e-(6i)e + (6a)e-(x) + 0-(x)0+(& 2 )©-(&i)) . 

Hence, 

PQ-(x) + Q-{x)P T - 0"(Px) 
= AQ-{x) + <d~{x)A T - Q~(Ax) 

- Q-(b 1 )Q-(x)Q-(b 1 ) - e-(b 1 )e-(x)e-(b 1 ) 

- Q-(b 1 )e+(x)Q-(b 2 ) - e-(6a)e+(x)0-(6i). (29) 
From ( I24db . it is thus evident that P obeys 

P0"(x) + Q-(x)P T - Q-(Px) = 0. 

Given that this equation is satisfied for all t, a direct 
application of Lemma [3] when t = shows that P = 0~(a) 
for a vector a as in (fTol l. which implies (12 lcb - 

Conversely, assume that J2 1 al >- (|2 lcl > are true. By Lemma 
|3] equations ( I24ab - (l24cb and ( I26cl i are automatically satisfied. 
Equation (I26db it is easily obtained by computing A + A T 
using d21cb - Furthermore, equations (I26at and (I26bt are a 
direct result of (ISfb . The validity of (|24db is shown by 
computing y40~(x) + _ (x)A T — _ (Ax). This calculation 
is analogous to the computation of d29l , but in this case it 
is already known that P = 0~(a). Thus, ( |24db is true given 
that d2 1 ab - (!2 lcb hold. To show that d26eb holds, one first 
computes Q + (Ax) using d21cb as follows 

Q + (Ax) = 6+ (0~(a)x) 

+ \ (0"(&i)0 + (e-(6i)x) - e+ (e-(6i)x) Q-{h) 
+e-(6 2 )e+ (e-(6a)x) - e+ (0-(& 2 )x) 0~(& 2 )) 

+ ^ (0-(& 2 )0 + (0 + (6i)x) - 0+ (0+(fti)x) e-(6a) 

-e-(6i)e+ (0 + (6 2 )x) + e+ (0+(& 2 )x) e-(&0) . 

Now, from the symmetry of + and ( [5g| , one observes that 

e+ (@ + (h)x) 

= e + (6 2 )6+(x) - 0-(x)0-(6 2 ) - -bfxI+-biX T 

n n 

= (e + {b 2 )e+{x) - 0-(x)©-(6 2 ) - -bfxi + -b tX T ) 

\ n n / 

= 0+(x)0+(6 2 ) - 0-(& 2 )0-(x) - -bfxl+-xbf. 

n n 

Therefore, 

e + (Ax) = 9+ (0~(a)x) 

+ 5 (0-(6i)0 + (x)0-(6i) - 0-(6!)0-(6i)0+(x) 

- 0+(x)0-(6!)0-(&i) + 0-(&x)0+(x)0-(6i) 

+ e-(& 2 )e+(x)e-(& 2 ) - 0-(& 2 )©-(& 2 )0 + (x) 
-e + (x)e-(fe 2 )e~(fe 2 ) + e-(fe 2 )e + (x)e-(6 2 )) 

+ \ f©-(6 2 )0+(6i)0+(x) - 0-(6 2 )0-(x)0-(6 1 ) 

- -e-(b 2 )blx + -Q-(b 2 )b 1 x T 
n n 

- e+(x)e + (b 1 )e~(b 2 ) + Q-(b 1 )e-(x)e-(b 2 ) 

+ -6^x0-(6 2 )--x6f©-(6 2 ) 
n n 

- e-(6 1 )e+(6 2 )6 + (x) + 0-(6i)0-(x)0-(6 2 ) 



+ -Q-(h x )blx - -e-(b 1 )b 2 x T 
n n 

+ Q+(x)e+(b 2 )Q-(b 1 ) - Q-(b 2 )Q-(x)Q-(b 1 ) 

- -bTxQ-(bi) + -xb^Q-(bi) \ . 
n n J 

Also, 

AQ+(x) + e + (x)A T 
= e-(a)Q+(x) + @+(x)Q-(a) T 

+ i (e-(&i)e-(6i)e+(x) + Q-(b 2 )Q-(b 2 )e+(x) 
+ e+(x)Q-(b 1 )e-(b 1 ) + e+(x)e-(b 2 )Q-(b 2 ) 
+ e-(fe 2 )e + (6i)e+(a;) - e + {x)& + {b 1 )e-(b 2 ) 

- e-(h)e + (b 2 )Q + {x) + Q + (x)e + (b 2 )e-(h)) . 

Since £ Q-{b 2 )b lX T = A x T and \ ^{b 2 )b x x T = xA%, It 
then follows that 

AQ~(x) + Q-{x)A T - Q+{Ax) 

= Q-(a)Q+(x) + @+(x)Q-(a) T - Q + (Q~(a)x) 

+ e-(6i)e+(x)e-(6x) + e-(b 2 )Q+(x)e~{b 2 ) 
+ @-(b 1 )Q-(x)Q-(b 2 ) + e-(6a)e-(x)e-(6i) 

- - (A x T + xAl) . (30) 
n 

Identity d5fb shows directly that 

0-(a)G+(x) + Q+(x)@-{a) T - e+(Q-(a)x) = 0. 
Hence, (f3Qb implies d26el ). which completes the proof. ■ 

Theorem 4: A physically realizable system satisfies the 
conditions of Theorem [3] 

Proof: Consider, without loss of generality, that n w = 1. 
Define C = \{C\ + iC 2 ). This means that C is playing 
the role of A in Theorem Q] By conditions ([m]) and (|mp of 
TheoremEl the fact that B x = 8~(6i) and B 2 = 6"(6 2 ) is 
automatic by considering b\ = C 2 and b 2 = C\. Next, one 
has from (0) of Theorem [2] that 

= - ( \{B X -iB 2 )(d +ic 2 y) = --ie-(c)cK 

n \2 J n 

Now, from d5d]i, Q~{Ao) is 

^e-(A)) = -e-(2ie-(c)c t ) 

= -2i (e-(C)6-(C t ) - 9" (tf)Q- (C)) . 
Since i?i is antisymmetric, it follows that 

BiBj - s 2 sf = -2i(e- (c)e _ ■(c t ) - e-(c t )e-(c)). 

Therefore, (I21bt holds. Finally, since the considered system 
is physically realizable then adding ( 1 12bb and its transpose 
gives 

= e+(6 a )e-(6 1 )-e-(6i)e + (6 2 ) 

+ e-(6a)e + (6i)-e + (6i)0-(6a) 

= |e+(A ), 



which concludes the proof. ■ 

The implication of this theorem is that a state space model, 
as given by dTOb and (Hit , describes an open ?i-level quantum 

system when it satisfies the physical realizability conditions 
since they also ensure preservation of commutation and 
anticommutation relations of SU(n). In addition, it is known 
that the physical realizability conditions can be employed to 
obtain the (S, L, H) parametrization of that quantum system. 

V. Conclusions 

Conditions for preserving commutations and anticommu- 
tation relations have been provided for QSDE's of the form 
( [Tol l. These results used explicitly the algebra of SU(n), 
and in particular algebraic expressions involving the anomaly 
tensor d had to be employed. Moreover, it was shown 
that the physical realizability implies the preservation of 
those relations, and therefore under physical realizability 
conditions the system given by ( TTOb and (fTTT > describes an 
open n-level quantum system. 
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